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Abstract 



This article is in continuation of our earlier article [37] in which computational solution of an unified 
reaction-diffusion equation of distributed order associated with Caputo derivatives as the time-derivative 
and Riesz-Feller derivative as space derivative is derived. In this article, we present computational solutions 
of distributed order fractional reaction-diffusion equations associated with Riemann-Liouville derivatives of 
p j ■ fractional orders as the time-derivatives and Riesz-Feller fractional derivatives as the space derivatives. The 

| method followed in deriving the solution is that of joint Laplace and Fourier transforms. The solution is 

derived in a closed and computational form in terms of the familiar Mittag-Leffler function. It provides an 
| elegant extension of the results given earlier by Chen et al. [1], Debnath [3], Saxena et al. [36], Haubold 

et al. [15] and Pagnini and Mainardi [30]. The results obtained are presented in the form of two theorems. 
Some interesting results associated with fractional Riesz derivatives are also derived as special cases of our 
findings. It will be seen that in case of distributed order fractional reaction-diffusion, the solution comes in 
a compact and closed form in terms of a generalization of the Kampe de Feriet hypergeometric series in two 
variables, defined by Srivastava and Daoust [46] (also see Appendix B). The convergence of the double series 
(*C) ' occurring in the solution is also given. 
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1. Introduction 



Distributed order sub-diffusion is discussed by Naber [26] . Distributed order fractional diffusion systems 
are studied, among others, by Saxton [40,41,42], Langlands [21], Sokolov et al. [43,44], Sokolov and Klafter 
[44], Saxena and Pagnini [38] and Nikolova and Boyadijiev [28], and recent monographs on the subject 
[4,20,22,25]. General models for reaction-diffusion systems are discussed by Wilhelmsson and Lazzaro [46], 
Henry and Wearne [17,18]. Henry et al. [19], Mainardi et al. [23,24], Haubold et al. [14,15,16], Saxena et 
al. [34,35,36,37,38] and others. Stability in reaction-diffusion systems and nonlinear oscillations have been 
discussed by Gafiychuk et al. [10,11]. Pattern formation in reaction-diffusion as well as non-Gaussian, non- 
Markovian, and non-Fickian phenomena related to astronomical, physical, chemical, and biological sciences 
can be found in the work of Cross and Hohenberg [2], Nicolis and Prigogine [27], and Haubold et al. [16]. 
Recently, Engler [6] discussed the speed of spread of fractional reaction-diffusion. 

In a recent article, Chen et al. [1] have derived the fundamental and numerical solution of a reaction- 
diffusion equation associated with the Riesz fractional derivative as the space derivative. Reaction-diffusion 
models associated with Riemann-Liouville fractional derivative as the time derivative and Riesz-Feller deriva- 
tive as the space derivative are recently discussed by Haubold et al. [14]. Such equations in case of Caputo 
fractional derivative are solved by Saxena et al. [37]. 

The main object of this article is to investigate the computational solutions of fraction reaction-diffusion 
equations (2.1) and (4.1). The results are obtained in a closed and computational forms. Due to general 
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character of the derived results, many known results given earlier by Chen et al. [1], Haubold ct al. [15] 
and Pagnini and Mainardi [30], Saxena et al. [37], and others, readily follow as special cases of our derived 
results. 

2. Solution of fractional reaction-diffusion equations 

In this section, we will investigate a computational solution of one-dimensional fractional reaction- 
diffusion equation (2.1) containing Riemann-Liouville derivative as the time-derivative and a finite number 
of Ricsz-Feller derivative as the space derivatives. The result obtained are in a compact and computational 
form in terms of the generalized Mittag-Leffler function, defined by (2.5) in the form of the following theorem: 

Theorem 1. Consider the following one- dimensional non-homogeneous unified fractional reaction- diffusion 
model associated with time-derivative D" defined by (Al), n g N, and Riesz-Feller space-derivatives 
x D%,..., x Dfc, defined by (A3): 



D°N(x,t)=^2^ x D%N(x,t) + <f>(x,t) (2.1) 

where t > 0, x g i?; a, 6i, 6 n , 71, 7„ are real parameters with the condition 

Hj > 0,0 < 7j- < 2,j = l,...,n, \6j\ < min (7,, 2 - 7,-), 1< a < 2; (2.2) 

l<j<n 

with initial conditions 

£>° _1 iV(a;,0) = f{x), D?- 2 N(x,0) = g(x), x g R, lim N(x, t) = 0, t > 0. (2.3) 

2:— >±oo 

Here oD^N(x, t) means the Riemann-Liouville fractional partial derivative of N(x, t) with respect to t of 
order v evaluated at t = 0,v = a — 1, a — 2; xDj 1 , X D2" are the Riesz-Feller space fractional derivatives 
with asymmetries 9i,...,9 n respectively. Further, o-Df is the Riemann-Liouville time-fractional derivative of 
order a; ...,\i n are arbitrary constants; f(x),g(x) and 4>{x,t) are given functions. Then for the solution 
of (2.1), subject to the above conditions, there holds the formula 

pa-l roo 



j-a-i roc 

N(x,t) = - / t*(k)E aa {-bt a )exp{~ikx)dk 

2lT J-oo 

/oc 
g*{k)E a . a _ 1 (-bt a ) exp{-ikx)dk 
-OO 



2vr 

1 /"* f°° 

+ 27 J * J ^*^^- E "<° cx P (-»fcx)dfc d& (2.4) 

where E a< p(z) is the generalized Mittag-Leffler function, defined by Wiman [48] (also see Erdelyi et al. [7], 
Dzherbashyan [5]) in the following form: 



z n 

E a .p{z) = ]T — Ty a,(3e C,R(a) > 0,K(/3) > (2.5) 

n=0 ^ n + 1 ' 



and 



n 



b=j2^<m- ( 2 - 6 ) 

Proof: In order to derive the solution of (2.1), we introduce the joint Laplace-Fourier transform int he form 

poo poo 

N*(k,s)= / c- st+lkx N{x,t)ddt 

JO J-00 
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where 3?(s) > 0, k > 0. If we apply the Laplace transform with respect to the time variable t, Fourier 
transform with respect to space variable x, use the initial conditions (2.2), (2.3) and the formula (A2), then 
the given equation transforms into the form 

n 

s a N*(k, s) - f*(k) - sg*(k) = - <J (k)N*(k, s) + 4>*(k, s) 

j=i 

where according to the conventions followed, the symbols ~ will stand for the Laplace transform with respect 
to time variable t and * represents the Fourier transform with respect to space variable x. Solving for N* (k, s) 
it yields 

Ar (M =m + ^ + ^ (2 . 7) 

s a + b s a + b s a + b 

where 

n 

&=I><2(*)- 

On taking the inverse Laplace transform of (2.7) by means of the formula 

= tP- a E^_ a+1 (-btP) (2.8) 



b + sP 

where SR(s) > 0, 5ft(/3) > 0, 5ft(/3 - a) > -1; |^| < 1 it is found that 

7V*(M) -* Q - 1 r(fc)^a,a(-W Q )+i Q -VW^a,a-l(-W Q ) 

+ / 0*(M-Or -1 £a,aHOd£- ( 2 ' 9 ) 

Jo 

The required solution (2.4) is now obtained by taking the inverse Fourier transform of (2.9). This completes 
the proof of Theorem 1. 

3. Special cases 

(i): If we take 6\ = ... = 6 n = then by virtue of the relation (All), we obtain the following corollary. 

Corollary 1.1. Consider the following one-dimensional non-homogeneous unified fractional reaction- diffusion^ 
model associated with time derivative o-D" obtained by (Al), n e N, and Riesz space derivatives jDjJ 1 , x-Do"! 
defined by (A3): 

n 

D?N(x,t) =^>j x D^N(x,t)+cf>(x,t), (3.1) 

where t > 0,i£ R; a, 71, 7„ are reo/ parameters with the constraints 

Hi > 0, <7j < 2,j = l,...,n,l < a < 2; (3.2) 

wii/j initial conditions 

D%- 1 N{x,Q) = f{x), D?-' 2 N(x,Q)=g(x), x&R, lim N{x,t) = 0, i > (3.3) 

x— >±oo 

w/iere i/ie various terms are as defined in (2.3). Then there holds the formula 

j.a—1 roo 

N(x,t) = / f*(k)E aa {-ct a )exp(-ikx)dk 

+ — — / g*(k)E a . a ^i(-ct a )cxp(-ikx)dk 
2tt J_ 00 

+ r _1 J <f>*(k,t-OE a , a (-ct a )eM-ikx)dkdli (3.4) 



where E a ^(z) is the generalized Mittag-Leffter function, defined by (2.5) and c = Y^j=i Mil^l 73 

Note that when g(x) = then Theorem 1 yields (2.9) where the middle term involving g*(k) will be absent. 

Corollary 1.2. Consider the same equation in (3.1) under the conditions 

Hi > 0, <7j < 2,j = l,...,n,\Bj\ < min (7^,2-7,-), < a < 1; (3.6) 

l<j<n 

with the initial conditions 

oD^N^x, 0) = 6(x), D?- 2 N(x, 0) = 0, x e R, lim N(x, t) = 0, t > (3.7) 

a;— ^±00 

w/iere 5(x) is Dirac de/ia function. Then for the fundamental solution there holds the formula 



N(x,t) 



t 



a-l r°° 



2tt 



/ 

J — C 



E a . a {—bt a ) exp(— ikx)dk 

(k, t - £)E ata (-b£ a ) exp(-ifcx)dfc d£ 



t pOO 

a-l 



(3- 



where E a> p{z) is the generalized Mittag-Leffler function defined by (2.5) and b is given in (2.6). 



If we further set <f) = 0, n = =11,61 = 6, 71 =7 and apply the formula (A14) then we have the 
following result: 



Corollary 1.3. Consider the following reaction- diffusion model 

D?N(x,t)=n x D2N(x,t) 
with the initial conditions 

oD^Nix, 0) = 6(x), x e -R, < a < 1, lim JV(z, t) = 0, i > 

cc— >±oo 

w/iere (1 is a diffusion constant, \i,t > 0;a,[3,6 are real parameters with the constraint 

< 7 < 2, \6\< min(7,2-7) 



(3.9) 
(3.10) 

(3.11) 



and S(x) is a Dirac delta function. Then the fundamental solution of (3.9) with initial conditions, there 
holds the formula 



N{x,t) = 



a-l 



7|x 



H-2,1 
■■"3,3 



(1,1/7), (a, a/7), (1,P)' 



(l,l/ 7 ),(l,l),(l,p) 



7 > 



(3.12) 



w/iere i?3' 3 (-) is the familiar H-function defined by (A12), also see [25]; P = \ 



2 a 



For 9 = 0, (3.12) reduces to the result given by Saxena et al. [38] in a slightly different form. On the 
other hand, if we further set a = 1 wc obtain the following result given by Chen et al. [1] in a different form: 



Corollary 1.4. Consider the following reaction- diffusion model 
d 



dt 

with the initial consitions 



N{x,t) = n x D%N(x,t) 



N(x, 0) = S(x), x £ R, lim N(x, t) = 0, t > 0, < (3 < 2, 

x— >±oo 



(3.13) 
(3.14) 



where n is a diffusion constant, fi,t > 0,(3 are real parameters and 5(x) is Dirac delta function. Then for 
the fundamental solution of (3.13) with initial conditions (3.14) there holds the formula 



N{x,t) 



tH. 



7R 



2.2 



\x\ 


(1, 1/7), (1, 1/2)- 




(1,1), (1,1/2) 



,7>0. 



(3.15) 



Remark 3.1. The result obtained by Chen et al. [1] is in terms of the Fourier transform of the Mittag- 
Lcffler function, whereas our result (3.15) is in terms of the H-function in a closed and computable form. It 
is interesting to observe that for n = 1, Theorem 1 reduces to one given by Haubold et al. [15]. 

4. Further results on distributed order reaction-diffusion systems 

In this section we will investigate a computational solution of a distributed order reaction-diffusion 
equation containing two Riemann-Liouville derivatives. The solution is obtained in terms of generalized 
Mittag-Lefner function due to Prabhakar [32]. The main result can be expressed in terms of Srivastava- 
Daoust hypergeometric function of two variables [45]. 

Theorem 2. Consider the following unified one- dimensional non-homogeneous reaction- diffusion equation 
of fractional order: 

n 

D?N(x,t)+a D?N(x,t) = x D%N(x, t) + (j>(x, t) (4.1) 

where t > 0, a, x € R; a, (3, 6\, 9 n are real parameters with the constraints 

Hj > 0, <7j < 2,j = l,...,n, \9 \ < min (7^,2-7^), 1 < a < 2, l</3<2 (4.2) 

l<j<n 

with the initial conditions 

D?- 1 N(x,0) = f 1 (x), a Dr 2 N(x,0)= 9l (x), o^f" 1 ^, 0) = f 2 (x), 

D^ 2 N{x, 0) = g 2 (x), xGR, lim N(x, t) = 0, t > (4.3) 

where the various quantities are as defined in Theorem 1, and f\{x), f 2 {x),g\{x),g 2 {x) and (fi(x,t) are given 
functions. Then for the solution of (4-1), subject to the above conditions (4-2), (4-3), there holds the formula 

j.a—1 poo 00 

N ^t) = - o- / (fi(k) +a/*(A : ))^(-a)'-^-«'-^+ 1 +(Q ^ )r (-&r)exp(- J fcx)dfc 

^~ J— 00 „_ n 



t 



r=0 

a — 2 roo 00 

/ (9*i(k) + ag* 2 {k)) Y.(- a y t(a ^ )r K^ +{a -f3)r-i(- hta ) eM~ikx)dk 

J—oo , — n 



2tt 



r=0 

roo 



i f-Z poo ^ poo 

+ / / ^-a)^"^ / <P%p,t- $E^ +{a _p )r (-bC)eM-ikx)dk d£ (4.4) 

^ J0 J-oo r= Q J — oo 



where ^(z) is the generalized Mittag-Leffler function ([32]) 



n=0 



where b = X)™ =1 Vj^ij (k). 

Proof: If we apply the Laplace transform with respect to the time variable t, Fourier transform with respect 
to space variable x and use the initial conditions (4. 2), (4. 3) and the formula (4.5), then the given equation 
transforms into the form 

s a N*(k,s) - f*(k) - sgl(k) + as^N*(k,s) - a/ 2 *(fc) - asg* 2 (k) 

n 

= -J2 ^% (k)N* (k, s) + 4>* (k, s). (4.6) 

3=1 

Solving for N*(k,s) it yields 



N*{k,s) 



mk)+afi{k) , S (gt(k) + ag*(k)) , 4>*(k) 



s Q + as/ 3 + 6 s a + as/ 3 + 6 s a + as? + b 



(4.7) 



where b is defined in (2.6). To invert the equation (4.7), we first invert the Laplace transform and then the 
Fourier transform. Thus to invert the Laplace transform we use the formula given in [40] : 



oP-l 



s a + as? + b 



t a - p T,(- a r t(a ^ )rE Z 1 H a -P)r- P+ i(- bta ) 



(4.8) 



r=0 



where SR(s) > 0, SR(a) > 0, SR(/3) > 0, 5ft(a - p) > -1, ft (a - /3) > 0, |^ 
of the Laplace transform to obtain 



< 1; and the convolution theorem 



Ar*(fc,i)=r- 1 [/r(fc) + a/|(fc)]^(-a)^-^+ 1 +(a _ /3)r (-6r) 

oo 

+ t«-*[ 9 l{k) + a«fe*(fc)] E(-°) rt(a ^ )r ^a 1 + (a- i8 )r-l(- 6 * a ) 
r=0 

oo - t 

+ 2(-«) r / ^(M-r^^CM^-^^ 

r=0 J ° 

Now, the application of the inverse Fourier transform gives the following solution: 



(4.9) 



N(x,t) 



j.a — 1 r°° 



2ir 



-OO 

±a — 2 poo 



/oo _ _ 

[f*(k) + ag* 2 {k)] £(-ar^-^i£+J +(a _^ r (-&n exp(-ifcs)dfc 



27T 



r=0 
oo 



[ 9l (fc) + a ff2 *(fc)]^(-a)^(«-«'-^+ 1 +(Q _ / , )r _ 1 (-&P)cxp(-*fc a; )dA: 



r=0 



/ \r ft rOC 

+ E^ L / r +(a - /3)r " 1 / 0*(M-£)expHfcr) 



K:: +{a -, )r (-bndk. 



Alternative form of the solution (4.4) 



(4.10) 



By virtue of the series representation of the generalized Mittag-Leffler function Ep (z) defined in (4.5), 
the expression 



r=0 



can be written as 



(l) r+u (-at a -r j y(-bt a ) u 



(4.11) 



OO OO 



r=0 u=0 



r\u\ T(a — p + 1 + (a — /3)r + au) 



j-a — p ol:0;0 

— 1 °1:0;0 



-at a -P-bt a 



[1:1;1]:-;- 

[a— p+l:a>— /3;a]: — ;- 



(4.12) 



where S'(-) is the Srivastava-Daoust generalization of the Kampe de Feriet hypergeometric series in two 
variables [46], for its definition, sec (B2) in Appendix B. Hence, Theorem 2 can now be stated in terms of 
the Srivastava-Daoust hypergeometric function of two variables in the following form: Under the conditions 
of Theorem 2, the unified one-dimensional fractional reaction-diffusion equation 
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D?N(x, t)+a D?N(x, t) = A x D}N(x, t) + <f>(x, t), 
has the solution given by 

t° 



(4.13) 



j.a-1 r°° 



N(x,t) 



2ir 
x <? 1:0;0 

x »~>1:0;0 



\fi(k) + afi(k)]exp(-ikx 

[1:1;1]:-:- 
[a-.ct — /3;a]: — : 



+ 



2tt 



at a -P,-bt a 

/OO 
-oo 



dfc 



ct-2 /-oo 



-at"-' 3 , -6t Q 



[a — l:a— /3;a] : 



dfc 



1 /"* Z" 00 

[i:i;i]-;- 

[a:a— /3;a]:- 

where SR(a) > 0, SR(/3) > 0, 3?(a - j3) > and b is defined in (2.6). 



x 9 1:0;0 

x J l:0;0 



dfc d£ 



(4.14) 



Note 4.1. By virtue of the lemma given in the Appendix B, the double infinite power series occurring in 
Theorem 2 converges for Sft(a) > 0, 3ft(a - /?) > 0. 

5. Special cases of Theorem 2 

If we set <7i(x) = 32 (x) = in Theorem 2 then it reduces to the following: 

Corollary 2.1.. Consider the following unified one- dimensional non-homogeneous reaction- diffusion equa- 
tion of fractional order: 



D?N{x, t) + a D^N{x, t) = ^ Mj xD%N(x, t) + cj>{x, t) 

w/iere i > 0, a, a; e ii; a, /3, #1, 0„, 71, 7„ are reaZ parameters with the constraint 

fij > 0, < 7j < 2, j = 1, n, 1^1 < min ( 7j , 2 - 7j ), 1< a < 2, 1< [3 < 2 

l<^<n 

wii/i t/ie initial conditions 

D?- 1 N(x,0) = f 1 (x), D?- 2 N(x,t) = 0, D^- 1 N(x,0) = f 2 (x), 
D^ 2 N{x,0) = 0, x e R, lim AT (a, i) = 0, t > 



(5.1) 



(5.2) 



(5.3) 



where the various quantities are as defined in Theorem 2. Then we have a special case of (4-10) and (4-12) 
with the corresponding changes. 

Other special cases of interest are the situations (i): n = 1, (ii): 0\ = ... = n = 0, (iii): fi(x) = f2{x) = 
5(x) in Theorem 2. We get the corresponding results from (4.10) and (4.14) by substitution. 
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Appendix A: Riemann-Liouville and Riesz-Feller fractional derivatives 

The Riemann-Liouville fractional derivative of order a > is defined as (Samko et al.[32, p. 37]; also 
see, Kilbas et al.[19]) 

1 d™ /"* 

D?f(x, t) = f(w _ q) — J (t - T ) n - a -\f{x, t)At, n = [a] + 1, n e N, t > (Al) 

where [a] means the integral part of the number a. 

The Laplace transform of the Riemann-Liouville fractional derivative is given by the following: (see 
Oldham and Spanier [29, (3.1.3)], Kilbas et al.[20], Podlubny [31], Samko et al.[33], Mainardi [22], Diethlm 
W) 

This derivative is useful in deriving the solutions of integral equations of fractional order governing certain 
physical problems of anomalous reaction and anomalous diffusion. In this connection, one can refer to the 
monograph by Dzherbashyan [5], Podlubny [31], Samko et al.[33], Oldham and Spanier [29], Kilbas et al.[20], 
Mainardi [22], Diethclm [4] and a recent paper on the subject [28]. Following Feller [8,9], it is convenient to 
define the Riesz-Feller space-fractional derivative of order a and skewness 9 in terms of its Fourier transform 
as 

F{ x D^f( X y,k} = -^ a (k)f*(k) (A3) 

where 

tP e a (k) = \k\ a exp[i(sign k)—], < a < 2, |0| < min{a, 2 - a}. (A4) 

Further, when 9 = 0, we have a symmetric operator with respect to x that can be interpreted as 

j2 i«/2 

(A5) 



n a — — 



dx 2 



This can be formally deduced by writing — (k) a = — (k 2 ) a / 2 . For < a < 2 and \6\ < minjo;, 2 — a}, the 
Riesz-Feller derivative can be shown to possess the following integral representation in x domain: 



c D%f(x) = £ii±^| sin[(a + 9)n/2] £ 



f(x + t)-f{x) 



l+a 



dxi 



+ sin 



m[(a-6)n/2] f 
Jo 



f(x-Q-f(x) 

j:l+a 



(A6) 
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For 9 = the Ricsz- Feller fractional derivative becomes the Riesz fractional derivative of order a for 1 < a < 2 
defined by analytic continuation in the whole range 0<a<2,a^l (see Gorenflo and Mainardi [12]) as 



where 



C D% = -\{I+ a - IZ a ] 
X = 1 - 



a r2-a 

dx 2 ± ■ 



2cos(mr/2)' ± 

The Weyl fractional integral operators are defined in the monograph by Samko et al. [33] as 



(A7) 
(AS) 



1 f°° 

(I + N)(x) = —Jj x - O^NiOdt (3 > 



(I*N)(x) = 



{i-xf^N{i)Ai, (3>0. 



(A9) 



Note Al. We note that x Dg is a pseudo differential operator. In particular, we have 



x Dl = — but x Dl * - 



For 9 = we have 



F{ x Dgf(x);k} = -\k\ a f*(k). 
The H- function is defined by means of a Mellin-Barnes type integral in the following manner [25]: 



(AW) 
(All) 



rjm,n ( \ rrm,n 

n <p,q \ Z ) p,q 



Tjm.n 



(b q ,B q ) 

(ai,A 1 ),...,(a p ,A p ) 



(6i,Bi),...,(6„B,) 



where i — v^T and 



e(0 



{ut=i r fe- + b jQ} {uu r ( x - a j_ M)} 



(A12) 



(413) 



and an empty product is interpreted as unity; m,n,p,q e N with < n < p, 1 < m < q, Ai, Bj £ R+, 
a,, bj € C, i = 1, -;P,j = 1, q such that the poles of T(bj + Bj£), j = 1, m are separated from those of 
T(l — cij — Aj£),j = 1, n, where N — {0, 1,2, ...}; R = (— oo, oo), R + = (0, oo) and C being the complex 
number field. A comprehensive account of the H-function is available from Mathai ct al. [25] and Kilbas et 
al. [20]. We also need the following result in the analysis that follows: Haubold et al. [14] has shown that 



F-^E^-at^k^x] 



a\x 



— H 3.3 



(atpy/ a 



(1,1/ 



a ),(7,0/a),(lp) 

(l,l/c0,(l,l),(l,p) 



(A1A) 



where 3?(a) > 0,3?(/3) > 0,3?( 7 ) > 0. 

Appendix B: Convergence of the double power series 

The following lemma given by the authors is needed in the analysis. 
Lemma. (Saxena et al. [39]). For all a,a,/3 > there holds the formula 

[1:1;1]: 



(l) m+n x m y n 



y y 

^o^o m!n! (°W^ 



t~i / \ c l:0;0 



x,y 



(Bl) 
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where S stands for the Srivastava-Daoust function (Srivastava et al. [4-5]), defined by (B2). 



Definition Bl. Srivastava-Daoust generalization of the Kampe de Feriet hypergeometric series in two 
variables is denned by the double hypergeometric series as [45]: 



nA:B;B' , x _ q *:a;a 
°C:D;D' V x i V) ~ °C:D;D 



_ qA:B;B' 



x,y 



EE 



[(a):0,#[(&):#[(6'):V>'] 

[(c):8,e]:[(d): V ]:l(d'y.r,'] 
B 



inii r( 0j + mOj + n^)}{nf =1 nb, + m^mnti m + nmx™y- 



rfn=o {IL=i r (cj + mSj 



™i)HIlf=i r(dj + ro»fe)HIlfIi IK + n V ')}m\n\ 



(B2) 



where the coefficients 9\, 6 a, ■ Vi, ^d' > 0- For the sake of brevity (a) is taken to denote the sequence 
of A parameters ai, a a with a similar interpretation for (b), (a"). Srivastava and Daoust have shown 
[45] that the series (B2) converges for all x,y e C when 



a = i + s j + E * - E °i - E *i > 



and 



a' = i+E £j + E^-E^-E^>o. 

j=i j=i j=i j=i 



(S3) 



(B4) 



For a detailed account of the convergence conditions of the double series, see Srivastava and Daoust [45]. 
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